Reversable heat flow through the carbon nanotube junctions 
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Microscopic mechanisms of externally controlled reversable heat flow through the carbon nanotube 
junctions (NJ) are studied theoretically. Our model suggests that the heat is transfered along the 
tube section T by electrons (e) and holes (h) moving ballistically in either in parallel or in opposite 
directions and accelerated by the bias source-drain voltage Vsd (Peltier effect). We compute the 
Seebeck coefficient a, electric a and thermal k conductivities and find that their magnitudes strongly 
depend on Vsd and Vg- The sign reversal of a versus the sign of Vg formerly observed experimentally 
is interpreted in this work in terms of so-called chiral tunneling phenomena (Klein paradox). 

PACS numbers: 73.23.Hk, 73.63.Kv, 73.40.Gk 



Physics of the heat transfer determines functionalit; 
precision and effectiveness of solid state nanocoolers[H_ 
y, [B[ which are environment friendly and have a lot of ap 
plications in the experimental physics, nanoelectronics, 
chemistry, industry and medicine. Therefore exploiting 
of new thermoelectric materials with high figures of merit 
Z ■ T (T being the temperature) attracts a lot of atten- 
tion. Recently such interest arose toward the carbon nan- 
otube and graphcnc junctions which electronic properties 
are highly unconventional [5| . The thermoelectric power 
experiments [I[ addressed single wall carbon nanotube 
junctions. In Ref. [f| a temperature difference AT in- 
duced a finite bias voltage AV P across the junction which 
sign changed versus the gate voltage Vg . A question here 
is how that unconventional thermoelectric behavior is re- 
lated to the intrinsic nature of the carbon nanotube and 
graphene? It is widely accepted that the charge carrier 
motion in carbon nanotubes and in graphene is essen- 
tially phase-correllated. For such a reason the conduct- 
ing electrons and holes in that materials behave as rela- 
tivistic massless 'chiral fermions' (CF) characterized by 
a 'pseudospin' (see review Q and references thereis). 

In this Letter we argue that the phase-correlated 
thermoelectric transport of charge carriers implicates a 
voltage-controlled and reversable heat flow through the 
single wall carbon nanotube junctions (Peltier effect). 
The enhancement of Z ■ T in those ID devices where 
the charge carriers propagate ballistically occurs due to 
strong van Hove singularities (VHS). The VHS position 
is tuned by the gate voltage Vq in respect to the Fermi 
level £f of the electrodes. When an VHS and £f match 
each other, it results in a sufficient density of charge car- 
riers which contribute to the electric conductivity despite 
the Fermi energy itself is relatively small. A finite gate 
voltage Vq ^ is not just merely supplies either the 
electrons or holes into T, but rather creates a poten- 
tial barrier (Uq > 0) or well (Z7 < 0) for chiral fermions 
transmitted across the junction as illustarted in Figs.[TJb- 
d). A finite sorce-drain d.c. bias voltage Vsd ^ applied 
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FIG. 1: (a) Typical carbon nanotube junction (NJ) controlled 
by the gate voltage Vg- The source (S) and drain (D) elec- 
trodes are attached to the nanotube section T, ii,2 are metal- 
lic stripes serving as termometers. (b) Schematics of the heat 
(Q) and electric (J) transport responsible for the reversable 
Peltier effect. (c,d) Energy diagrams of basic microscopic tun- 
neling processes involving electrons (e) and holes (h) across 
the chiral barrier (c) and well (d). 



across the junction tilts the chiral barrier and causes the 
electric charge carriers to accelerate. This means that 
the energy of the charge carriers changes by 5 — eVsD, 
which inflicts a local temperature change AT at the cor- 
responding junction's end. The sign of 6 (and hence of 
AT) depends on whether the charge carriers inside the 
barrier are electrons or holes. Other thermal character- 
istics of the carbon nanotube and graphene junctions de- 
pend on precise microscopic mechanisms acting on the 
nanoscale. The setup shown in Fig. Q] allows controlling 
the sign of AT by merely flipping the sign of the gate 
voltage [see the sketch in Fig. [Tfb)]. The source (S) and 
drain (D) electrodes can either be metallic or made of 
the nanotube/graphene itself. The potential chiral bar- 
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rier Uq is induced by the gate voltage Vq from a Si gate. If 
Vq > 0, the transport mechanism is the chiral tunneling 
(CT). It presumes a constructive quantum interference 
between an incoming electron (e) and a hole (h) mov- 
ing inside the graphene barrier in a reverse direction and 
is characterized by the same pseudospin as sketched in 
Fig. [He) (Klein paradox^, Q). The interference pattern 
resulting from the chiral tunneling is very sensitive to the 
phase difference ip between the e and h wavefunctions. 
During the CT process the electric current is transferred 
by holes, which ballistically propagate in a reverse di- 
rection. That means the heat and the electric current 
are directed in antiparallel. When the gate voltage is re- 
versed (Vg < 0), the heat and the electric current flow in 
parallel, since both of them are transferred by electrons. 
This reverses the heat flow along T as compared to the 
former case Vq > 0. We emphasize that the reversable 
heat flow originates from the phase correlation between 
the electrons and holes which constitutes its quantum 
mechanical origin. Our finding suggest that the thermo- 
electric transport serves as another independent probe of 
the phase-coherence in the carbon nanotube junctions. 

Thermoelectric transport properties of the Peltier cool- 
ers are characterized by the figure of merit Z ■ T. For a 
symmetric setup shown in Fig. [2lc) Z takes the form 
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In Eq. the lattice thermal conductivity is kl = 

ClVlIl/^, Gl is the graphene heat capacity, which at 
T = 300K is C L = 8.5 J/(mol-K), v L w O.lv is the 
graphene sound velocity, where v = 8.1 • I0 5 m/s is 
the Fermi velocity in graphene, II is the phonon free 
path inside the nanotube. A most conservative esti- 
mation used here implies that II is comparable with 
the tube diameter, i.e., Il ~ d^. The Seebeck coeffi- 
cient of electron (+) and holes (-), a±, is obtained as 
a± = —l/(q±T)(L^/L^) , where q± = =Fe is the elec- 
tric charge of the electron (hole) . Besides in Eq. (JTJ (7 scq 
is the conductivity of two equal size conducting tubes T 
connected in a sequence, cr scq = o + <j-/ (a + + a J) where 

a± = L± and K njP is the electron (hole) thermal conduc- 
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/l4 0) ) /(e 2 T). The func- 



tivity k± = yL 
tions L± entering the above formulas are 



-(«) 



2e 2 iV(0) / (e-V SD ) a \tf\ 2 {-^)ds (2) 



The functions are computed from a microscopic 

model. If one neglects by detailed energy dependence 
of tf and merely sets \tf\ 2 -> T, then a p . n = e 2 N(0)T 
where N (0) is the electron density of states at the Fermi 
level of electrodes, and 
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a± 



Tlog2-F SD /2 



q ± T tt 2 T 2 + W sd (Vsb - 4T log 2) 



(3) 



and the electron (hole) thermal conductivity n n _ p be- 
comes 

k± = N (0) T(7r 2 T 2 /3 + V SD (V SB - 4Tlog2) 

-4(Tlog2-VW2) 2 )/T (4) 

The estimations using Eqs. (HJ)- (J3J) give the figure of 
merit for a symmetric setup as 



T ■ Z 



288 (log2-y SD /2) 2 
ot 2 z Pz + e 2 T • kl 



(5) 



where a z = n 2 T 2 + 3^ sd (^sd - 4Tlog2), Z = (tt 2 T 2 + 
3^sd(^sd - 4Tlog2))/3 - 4(Tlog2 - V SB /2) 2 . One 
may notice that T ■ Z does not depend on the prop- 
erties of electrodes but instead strongly depends on T 
and Vsd- For kl = 10 the product T ■ Z may be in- 
creased by 2 orders of magnitude when neglecting by the 
phonon heat conductivity. In Eq. ([21)-© we used that 
L^ 0) = e 2 iV(0)T L ( l ] = 2e 2 N{Q)f ■ (Tlog2 - Vsd/2), 

and L { ±> = e 2 N(0)f(ir 2 T 2 /3 + V S u(V SB -4Tlog2)). 
The above rough estimations illustrate the temperature 
behaviour of the figure of merit disregarding however the 
microscopic mechanisms of the thermoelectric transport. 
A deeper understanding is achieved with using of a mi- 
croscopic model where the electron and hole envelope 
wave functions satisfy the Dirac equation. That model 
allows computing the transparency T = \tf\ 2 of the chi- 
ral barrier and well shown in Fig. HKb-d). We calculate 
the steady state thermoelectric characteristics of an NJ 
with metallic electrodes controlled by the gate voltage 
Vq and biased by the source-drain voltage Vsd- ln par- 
ticular we will see that the gate voltage reverses the heat 
flow along T, which is well consistent with the earlier 
experiments [5] . The figure of merit Z ■ T strongly in- 
creases at certain magnitudes of Vg and Vsd ■ The source 
drain bias voltage Vsd in our setup drops entirely on the 
nanotube [see its profile below Fig. 1(b)], which corre- 
sponds to a trapezoidal shape of the chiral barrier tilted 
proportionally to Vsd- For the NJ having finite dimen- 
sions, the motion of CF is quantized. The quantization 
imposes additional constrains on the directional tunnel- 
ing diagram. The Dirac equation for fermions is written 
as 

- ihv ((a x ®l)d x + (a y ® f g ) d y ) * + eU (x) * = 

where v = c/300 is the massless fermion speed, e is the 
CF energy, &i and fk are the Pauli matrices, {i,k} — 
1 ... 3, the barrier potential U(x) is induced by the gate 
voltage Vg- In the steady state, when the EF is off, 
the electric current is fully suppressed when Vsd < Uq 
(for typical gate voltage Vg = 1 V and the Si02 thick- 
ness d = 300 nm one finds [l(| Uq — 2 meV). The 
CF energy in the leads (x < and x > L) reads 

e = eVoo ± hv^k 2 + q 2 (n) = Hv (k Fm ± ^Jk 2 + q 2 (n)^j 
where kp m is the Fermi wavevector in the S and D 
electrodes. Then in the electrodes one takes k = 
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y/ (e/hv — kFm) 2 — 1v {n) 2 ■ Inside the tube one uses 

k = J ((e — eU )/hv) 2 — q v (n) 2 , where the electron 
wave vector in the transverasal direction is q u (n) — 
(2 /<1t) {n — f/3), is the nanotube diameter, n is the 
electron subband index, v = n + m — 3N, n, m, and 
N are integer numbers related to the translation vector 
C = na + mb of the graphene lattice, a and b are prim- 
itive translation vectors. In particular index v = when 
the [to, n] tube is metallic while v ^ for semiconducting 
and dielectric tubes. 

In a simplest case Vsd = from continuous boundary 
conditions one gets 

t n = 2e~ iip kks~s/V (7) 

where T> = 2kkss cos(kL) — iOsin(fcL), L is the tube 

length, K = \jk 2 + g 2 , 8 = kk(s 2 + s 2 ) — 2ssq 2 , k — 

y 'k 2 + q 2 , s = 1, s = sign (e — eUo). The transmission 
amplitude t e is computed using a simple approach, which 
represents the CF envelope wavefunctions ^ as plane 
waves for a rectangular barrier and via Airy functions [l5j 
for a trapezoidal barrier. In particular, the scattering 
state for a trapezoidal barrier is constructed as 

* = 6»(-x)[xn,fcAi(fc, x) + r n Xn-kBi(k, x)] 
+9(x - L)i„x„, fc Ai(fc, x - L) + 0(x)0(L - x) 
x [a nXnt ~ k Ai(k, x)e* + /3„x„,-fcBi(fc, x)e^} (8) 

where Ai(fc,a:) and T5i(k,x) are the Airy functions, 
and we also introduced the auxiliary functions 
Xn,k(y) = flnlT) ® (\l)+z n .k\l))e^y +a' n \i) (8) 
(^,fc|T) + l4»e l9 " y +i«IT> ® (2^ jfc |T> + U» e-^" 
+ b n\i) ® (IT) +2„, fc ||))e- l< ?"« where (gi means the 
Kronecker product, ||) T = ( 1 ) and \l) T = ( 1 ) 
are 1x2 matrices, z n ^ — ±(k + iq n ) / \Jk 2 + q 2 , where 
± signs apply to conductive (valence) bands, k is positive 
for conductance band and negative for the valence band, 
the factor z n ^ satisfies the identity z nt kZ nt -k = — 1. 
Then one implements the above Eqs. dT|)- ([2|) , and {J} to 
obtain a±, a±, «±, and T • Z. 

Our theoretical model interpretes the experimental 
data of Ref. Q in terms of chiral tunneling. We find that 
the electric and heat transport characteristics of the nan- 
otube junctions are pretty much determined by details of 
their electron spectrum. In particular, the electric and 
heat transport characteristics of the semiconducting nan- 
otube (with v = 1) junctions shown in Fig. [2] for two 
SWCNT diameters (Lt = 0.75 (dashed curves) and c2t = 
1.5 (solid curves) computed at T — 0.5 remarkably de- 
pend on the barrier height Uq. The d.c. electric conduc- 
tivity a+(Uo) plotted in units of e 2 v F N (0) in Fig. [2^a) 
shows a huge dip in vicinity of U$ ~ 0, which comes 

from the dispersion law k oc \J (eUo/hv) 2 — (2z//3<ir) 2 
in the lowest subband with n — 0. The upper electron 
subbands with n = 1,2,3 give a fine structure in the 




FIG. 2: (a) The d.c. electric conductivity <r+(f/o) [in units of 
e 2 v F N (0)] for the two SWCNT diameters d T = 0.75 (dashed 
curve) and dx = 1-5 (solid curve) computed at T = 0.5 and 
v = 1. (b) The net Seebeck coefficient «+ + a_ [in units of 
l/(k B TeVso)] versus Uo for the same SWCNTs at T = 0.5 
and v = 1. Indices 0, 1, 2 in (b) mark features coming from 
the electron energy subbands with n — 0, 1, 2 correspondingly, 
(c) Schematics of the reversable Peltier nanocooler combined 
of two NJs [see the former Fig. 1(a)] with individual gate volt- 
age Vq 1,2 ' control each, (d) The corresponding dependence 
for the figure of merit T ■ Z(Uq). 



transport characteristics which is remakably pronounced 
in all the a±(Uo), (x±(Uq), k±(Uq), and T- Z(Z7 ) curves, 
as is evident from Fig.O To distinguish the contributions 
from different electron subbands with n — 0, 1, 2 we mark 
the corresponding features by indices 0,1,2 in Fig. [2b). 
The positions and magnitudes of those features reflect 
the electron band structure of the tube or the graphene 
stripe. Similar fine structure features coming from differ- 
ent electron subbands are visible also in Figs. [2a), (c,d) 
for two the SWCNT diameters (It = 0.75 (dashed curve) 
and c?t = 1.5 (solid curve) computed at T = 0.5 and 
v = 1. The reversable heat flow due to the chiral tunnel- 
ing is well pronounced in the Seebeck coefficient a±. In 
Fig. [2(b) we plot the net Seebeck coefficient a + + a_ [in 
units of 1/ (k&TeVsT>)]- One see that a+ + a_ reverses 
its sign as the sign of Uq changes. This sign change along 
with the specific one-dimensional bandstructure of the 
single wall carbon nanotubes can be exploited for creat- 
ing of very efficient Peltier nanocoolers. An example of 
that cooler setup is sketched in Fig. [2jc) . The reversable 
Peltier nanocooler is combined of two NJs [see the for- 
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mer Fig. 1 (a)] with individual gate voltage ' control 
for each of the tube. A significant enhancement of the 
figure of merit Z ■ T in those ID devices is achieved by 
matching of the van Hove singularities (VHS) with posi- 
tion of the Fermi level £f in the electrodes by changing 
the gate voltage Vq and the bias voltage Vsd ■ The afore- 
mentioned peculiarities coming from the chiral tunneling 
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FIG. 3: The electric conductivity a±, Seebeck coefficient a±, 
the thermal conductivity k± , and the figure of merit Z of the 
reversable Peltier nanocooler at T = 0.5 versus the source- 
drain voltage. All the units are the same as in former Fig. 
2. 

and from the electron band structure actually determine 
behavior of T-Z versus Uq shown in Fig. Eld). From that 
Fig- [Hid) one can see that for the semiconducting [y = 1) 
tube with diameter d = 0.75 the product T ■ Z may ex- 
ceed 1. Further improvement of the T ■ Z product about 
by one order of magnitude is acomplished by applying 
a finite bias voltage Vsd as shown in plots of cr±(Vsr>), 
a±(Vsn), k±(Vsd), and T ■ Z(Vsd) in Fig. EM). The 
transport coefficients were computed at the temperature 
T = 0.3 (in units C/ , which for Uq = 50 mV corresponds 
T = 150 K). One may notice that the T ■ Z product also 
depends on the tube diameter dj>. For a thinner tube 
with d = 0.75 one gets T ■ Z = 7 already at Vsd = 1.2 
while for the thicker tube d = 1.5 one achieves T ■ Z = 8 



by applying a much higher bias voltage Vsd = 6.5. Simi- 
lar double peak structure T • Z(Vsd) coming from match- 
ing of the VHSs with ep of electrodes [see solid curve in 
Fig. Eld)] had also been proclaimed in Bi nanowires [1] . 
More significant increase of the T ■ Z product is obtained 
at lower temperatutes. Using the cooler setup shown in 
Fig. Efc) one may potentially get even T ■ Z = 10 2 — 10 3 
for temperatures T = 10K down to 1 K and with an 
appropriate tuning by Vg and Vsd simultaneously. In- 
clusion of the Shottky barriers and Coulomb blockade 
in the present model is pretty straiforward and will be 
given elsewhere. Our study indicates that narrow semi- 
conducting single wall CNTs are more efficient in coolers 
than metallic CNT and graphene stripes (where there is 
no VHS for the lower subband n = 0). An evident rea- 
son is that the matching of VHSs with ep in electrodes 
results in a better cooling performance. Another strong 
benefit of the NJ cooler is that one avoids a mutual com- 
pensation of a+ and cv_ in Eq. [T]as it happens, e.g., for 
Bi nanowire cooler (see Ref. [4j) since one appropriately 
flips signs of a± merely by changing Ve- 
in conclusion we suggested a theoretical model of the 
electric and heat transport across the carbon nanotube 
junction. We find that the heat transport strongly de- 
pends on the polarity and magnitude of the source-drain 
and gate voltages. That voltages tune the basic trans- 
port properties of the carbon nanotubes via affecting of 
their electronic structure and phase-correlated transport 
of charge carriers. Due to resonant character of chiral 
tunneling and low inelastic scattering rates the heat cur- 
rent density is typically much higher than in ordinary 
semiconducting devices. At the same time, the Peltier 
effect is well pronounced up to relatively high tempera- 
tures. The reversable Peltier effect shows great promises 
for cooling nanodevices and termometers. 
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